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Abstract. We conduct several sets of simulations with a ver- 
sion of NASA’s Goddard Earth Observing System, version 5, 
(GEOS-5) Atmospheric Global Climate Model (AGCM) 
equipped with a two-moment cloud microphysical scheme to 
understand the role of biomass burning aerosol (BBA) emis- 
sions in Southeast Asia (SEA) in the pre-monsoon period of 
February-May. Our experiments are designed so that both di- 
rect and indirect aerosol effects can be evaluated. For clima- 
tologically prescribed monthly sea surface temperatures, we 
conduct sets of model integrations with and without biomass 
burning emissions in the area of peak burning activity, and 
with direct aerosol radiative effects either active or inac- 
tive. Taking appropriate differences between AGCM exper- 
iment sets, we find that BBA affects liquid clouds in sta- 
tistically significantly ways, increasing cloud droplet num- 
ber concentrations, decreasing droplet effective radii (i.e., a 
classic aerosol indirect effect), and locally suppressing pre- 
cipitation due to a deceleration of the autoconversion pro- 
cess, with the latter effect apparently also leading to cloud 
condensate increases. Geographical re-arrangements of pre- 
cipitation patterns, with precipitation increases downwind of 
aerosol sources are also seen, most likely because of advec- 
tion of weakly precipitating cloud fields. Somewhat unex- 
pectedly, the change in cloud radiative effect (cloud forcing) 
at surface is in the direction of lesser cooling because of de- 
creases in cloud fraction. Overall, however, because of di- 
rect radiative effect contributions, aerosols exert a net nega- 
tive forcing at both the top of the atmosphere and, perhaps 
most importantly, the surface, where decreased evaporation 
triggers feedbacks that further reduce precipitation. Invok- 
ing the approximation that direct and indirect aerosol ef- 


fects are additive, we estimate that the overall precipitation 
reduction is about 40 % due to the direct effects of absorb- 
ing aerosols, which stabilize the atmosphere and reduce sur- 
face latent heat fluxes via cooler land surface temperatures. 
Further refinements of our two-moment cloud microphysics 
scheme are needed for a more complete examination of the 
role of aerosol-convection interactions in the seasonal devel- 
opment of the SEA monsoon. 


1 Introduction 

Use of fossil fuels for ever-growing energy demands, particu- 
larly in developing countries, has led to increased concentra- 
tions of aerosol-laden combustion by-products, especially in 
the planetary boundary layer (PBL) (Roelofs, 2013). Moor- 
thy et al. (2013) estimate that aerosols over India have been 
increasing at the rate of 2-4 % per year over the last three 
decades, resulting in doubled aerosol optical depth (AOD) in 
India’s lower atmosphere. Similar changes are expected over 
other regions such as Southeast Asia (SEA). Biomass burn- 
ing (BB) is an age-old method of disposing of agricultural 
trash (Taylor, 2010) and in SEA, it occurs primarily during 
the spring season (i.e., February-March- April, FMA; Gau- 
tam et al., 2013). Over SEA, the combustion by-products re- 
leased into the atmosphere contain large quantities of bio- 
genic aerosol/carbon particles whose quantitative estimates 
are being tabulated with extensive measurements (Wiedin- 
myer et al., 2011). 

Biomass burning aerosol (BBA) can affect the atmo- 
spheric circulation in several ways. BBA absorbs and reflects 
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Table 1 . Experimental designs for our GEOS-5 AGCM simulations. 
“Zero” stands for zero BB emissions over the green dash box region 
of Fig. la, “High” for high (year 2007, Fig. lb) emissions. Symbols 
under aerosol direct effect and indirect effect indicate experiments 
with (O) and without (X) the effect. 



BB 

Direct effect 

Indirect effect 

HighBoth 

High 

O 

O 

ZeroBoth 

Zero 

O 

O 

Highlnd 

High 

X 

o 

Zerolnd 

Zero 

X 

o 


solar radiation, thereby reducing the solar radiation reaching 
the surface, reducing surface sensible and latent heat fluxes 
(Ramanathan et al., 2005). Chung and Ramanathan (2006) 
showed the so-called “Atmospheric Brown Cloud” decreases 
the surface solar radiation flux, which reduces surface evap- 
oration while also weakening latitudinal sea surface temper- 
ature (SST) gradients and stabilizing the troposphere, caus- 
ing monsoon rainfall decreases. On the other hand, absorp- 
tion of solar radiation at the aerosol level warms the local 
atmosphere, inducing elevated heating that can invigorate 
air mass convergence near the surface and, with the addi- 
tion of sensible and latent heat, can make the PBL unstable 
enough to promote moist convection (Lau et al., 2006; Lau 
and Kim, 2006, 2014). The net outcome of the resulting com- 
plex feedback interactions may either increase or decrease 
local rainfall (Meehl et al., 2008). Furthermore, many par- 
ticles from BB emission are active cloud condensation nu- 
clei (CCN) (Petters et al., 2009). Hence more BB emission 
leads to more CCN and ice nuclei (IN) and thereby more 
cloud particles. If we assume that the net condensate pro- 
duction is solely governed by cloud-scale dynamics, more 
CCN would imply larger numbers but smaller cloud droplets 
and thereby an increase in cloud albedo (Twomey, 1977). 
Smaller cloud particles would also hamper the autoconver- 
sion of cloud water into precipitation, so the presence of BB 
sources is expected to reduce the precipitation production 
rate and increase cloud lifetime (Albrecht, 1989). However, 
this process is only applicable to warm rain. Observations 
show that cold and mixed cloud regimes have complicated 
responses, as summarized in Tao et al. (2012, Table 1). Li et 
al. (2011), with extensive observational analysis in the At- 
mospheric Radiation Measurement (ARM) site at Southern 
Great Plains, show that cloud-top height and thickness in- 
crease with aerosol concentration in mixed-phase clouds and 
rain increases with aerosol concentration in deep clouds, but 
declines in clouds that have low liquid-water content. 

Satellite data reveal that during FMA, the SEA region ex- 
hibits the highest aerosol concentrations, an order of mag- 
nitude greater than that in the summer monsoon May-June- 
July -August (MJJA) season because of more BBA sources 
during dry FMA (Ichoku et al., 2008; Lin et al., 2009) and 


less wet scavenging of aerosols compared to rainy MJJA. Ac- 
cordingly, aerosol optical depth and aerosol-activated cloud 
particle numbers are expected to be much larger in FMA than 
MJJA. This is the main reason for focusing this investigation 
of BBA direct and indirect effects on FMA and the transition 
month of May. Our working hypothesis is that high aerosol 
number concentration in FMA has a strong influence on the 
radiative forcing, circulation, and precipitation of the local 
and surrounding region. Important factors woidd be aerosol- 
cloud-radiation interactions. The cloud cover over SEA is 
generally composed of stratiform, low-altitude clouds asso- 
ciated with frontal systems that originate in China (Hsu et 
al., 2003). The major type of precipitation would be warm 
rainfall in the focused season. 

Direct and indirect effects of aerosols are intrinsically in- 
teractive, and therefore their combined effects can be very 
different from their linear sum. Even though the fundamen- 
tal physics of aerosol direct and indirect effects is reasonably 
well understood, uncertainty of aerosol data under cloudy 
conditions and complexities in coupling the aerosol-cloud- 
radiation interactions prohibit a better understanding of the 
impact of these processes. For example, a positive relation- 
ship between AOD and total cloud cover (TCC) was shown 
in satellite data (Kaufman et al., 2005; Kaufman and Koren, 
2006), but the dominant contribution to the AOD-TCC re- 
lationship have been attributed to aerosol swelling in humid 
air rather than the direct effects of aerosols on the cloud fields 
(Quass et al., 2010). 

Aerosol-cloud interaction effects on South Asia to East 
Asia circulation and monsoons have been the subject of in- 
vestigations with regional models (Wu et al., 2013; Lim et 
al., 2014) as well as global climate models (Bollasina et 
al., 2011; Ganguly et al., 2012; Guo et al., 2013). Our cur- 
rent cloud physics scheme, Microphysics of clouds with Re- 
laxed Arakawa-Schubert moist convection upgraded with 
prognostic aerosol-cloud interactions (McRAS-AC; Sud et 
al., 2013), has indirect effect simulation capabilities, and has 
been implemented in the GEOS-5 AGCM (Rienecker et al., 
2008). It provides an opportunity to perform simulation stud- 
ies to systematically assess the influence of BBA on rainfall 
and circulation in SEA. Clearly, constrained model simula- 
tions are one plausible way to better distinguish between the 
roles of direct and indirect effects and their interactive influ- 
ences that depend on circulation, cloud types, and aerosol- 
dependent cloud microphysics. While in principle these ef- 
fects can be properly simulated only with a coupled ocean- 
atmosphere model, as a first step we use an AGCM with pre- 
scribed monthly SSTs and with aerosol emission anomalies 
prescribed from the Quick Fire Emissions Dataset (QFED) 
data set (see Sect. 2. 1). This way we can isolate the influence 
of BBA over land by comparative assessments of circulation 
and rainfall changes in neighboring regions. 

In this endeavor, we perform a comprehensive model sim- 
ulation study with the physically interactive aerosol-cloud- 
radiation treatment of McRAS-AC as implemented in the 
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GEOS-5 AGCM, in order to better understand the spatiotem- 
poral modulation of the SEA pre-monsoon season by BBA. 
Section 2 describes the data set, model and experimental de- 
sign; results and a summary are presented in Sects. 3 and 4, 
respectively. 

2 Data, model and experimental design 

2.1 Data sets for aerosol effect analysis 

BB is a major source of primary emissions of carbonaceous 
aerosols over the SEA region. QFED (Darmenov and da 
Silva, 2013) was developed to meet the needs of the NASA 
Goddard Earth Observing System Model (GEOS) with re- 
gard to atmospheric constituent modeling and data assimila- 
tion of BB events. QFED is based on global gridded fire ra- 
diative power, derived from the Moderate Resolution Imag- 
ing Spectroradiometer (MODIS) Level 2 fire product. QFED 
is used not only as a BB inventory for the global Goddard 
Chemistry Aerosol Radiation and Transport (GOCART, Chin 
et al., 2002, Colarco et al., 2010) model in the GEOS-5 sys- 
tem, but also as an index indicating high BB days for our 
composite analysis. Version 2.2 used in this study covers the 
period from January 2003 to December 2010. The QFED 
Level-3 products are available at 0.3125 x 0.25°horizontal 
resolution, but are degraded to 2.5 x 2.0 degree for use in the 
present model simulation. The 1° MODIS Aqua level 3 daily 
product (MYD03 D3) is used for aerosol optical depth (Chu 
et al., 2002) and liquid cloud effective radius (Platnick et al., 
2003). The data cover the period from July 2002 to present. 
For precipitation, 1° daily Global Precipitation Climatology 
Project (GPCP-1DD; Huffman et al., 2001) data are used, 
covering the period October 1996 to present. 

2.2 GEOS-5 AGCM with double moment microphysics 
and updated radiation 

The numerical model used for this study is the GEOS- 
5 AGCM, version Fortuna 2.5 documented by Molod 
et al. (2012). In the current application, McRAS-AC re- 
places the cloud scheme of the baseline model. McRAS- 
AC synthesizes the initial version of McRAS (described in 
Sud and Walker, 1999, 2003) with subsequently developed 
aerosol-cloud interaction microphysics described in Sud and 
Lee (2007). The latest modification to McRAS-AC includes 
the addition of Barahona and Nenes (2009a, b) ice nucle- 
ation for mixed phase and ice phase clouds, as well as Foun- 
toukis and Nenes (2005) liquid droplet formation. The pre- 
cipitation parameterization remains as before, namely Sud 
and Lee (2007) for the liquid phase and Sundqvist (1988) 
for the mixed and ice phases. In-cloud evaporation, precip- 
itation and self-collection of cloud water are parameterized 
according to Sud and Lee (2007), employing a reformu- 
lated version of the Seifert and Beheng (2001, 2006) param- 
eterization to handle the much thicker cloud layers encoun- 


tered in a coarse resolution GCM. These algorithms work 
seamlessly across widely varying vertical model-layer thick- 
nesses. Any change in the cloud water substance mass by 
condensation/deposition and/or collection by precipitation 
works interactively through an implicit backward numerical 
integration that approximates the solution of the basic non- 
linear coupled differential equations for the cloud source and 
sink terms of the mass balance tendency equation. Despite 
using the observationally based Sundqvist (1988) equations 
for the mixed phase and ice phase precipitation tendencies, 
the implementation of the Barahona and Nenes (2009a, b) 
ice nucleation and the Bergeron-Findeisen cloud water-to- 
ice mass transfer (Rotstayn et al., 2000) allows for a reason- 
able separation of cloud liquid and ice mass fractions with 
their accompanying liquid and ice particle number concen- 
trations. Homogenous freezing of in-cloud liquid droplets 
surviving below — 38 °C is enforced by assuming instanta- 
neous freezing. Aerosol-cloud interactions are implemented 
into both stratiform (large-scale) clouds, and convective tow- 
ers topped by detraining convective anvils that transform 
into large-scale clouds at a prescribed timescale of an hour. 
Sud et al. (2013) provides a much more comprehensive dis- 
cussion of McRAS-AC (including treatment of the differ- 
ent cloud types) and its comparative performance against 
the cloud scheme of the baseline model. That study also in- 
cludes sensitivity studies with an interactive aerosol module 
and modified aerosol size distribution. The model used in 
the current study contains all the upgrades outlined above. 
The CFMIP Observation Simulator Package (COSP, http: 
//cfmip.metoffice.com/COSP.html) is also employed online 
in our experiments. Because of the significant dilferences 
between the way clouds are observed and represented in 
models, a “satellite simulator” facilitates proper compari- 
son and validations of the key simulated cloud and radia- 
tion fields against observations (Klein et al., 2013). The GO- 
CART module provides prognostic aerosols fields consist- 
ing of five aerosol species with fifteen modes. There are five 
modes of dust and sea salt sorted in dilferent particle size 
bins; there are two modes of organic and black carbon to sort 
hydrophilic and hydrophobic particles; and one mode of sul- 
fate particles. All aerosol modes are assumed to be “exter- 
nally” mixed. The GOCART module runs interactively and 
provides prognostic aerosol fields within the AGCM. 

Accurate radiation calculations are also very important for 
properly simulating aerosol direct/indirect effect. Our means 
of calculating realistic cloud radiative effect (CRE) is the ad- 
vanced RRTMG radiative transfer package (Clough et al., 
2005) equipped with a subcolumn generator in the GEOS- 
5 AGCM. RRTMG can be run in Monte Carlo Indepen- 
dent Column Approximation (McICA) mode (Pincus et al., 
2003) that operates on subcolumns with either clear or com- 
pletely overcast cloud layers produced by a cloud generator. 
Whether the cloud condensate in a particular layer is dif- 
ferent from subcolumn to subcolumn depends on the spe- 
cific assumptions about horizontal cloud heterogeneity as 
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Figure 1 . QFED BB emission data (OC data is used for this figure, unit of ggm -2 s _1 ). (a) February-March-April mean and (b) averaged 
time series for dashed box area. Time series are smoothed by a 7-day moving average. Orange line represents data for year 2007 for high BB 
simulation experiments. Other years are plotted in gray. 


determined by distributions of condensate specified within 
the cloud generator. A prior implementation of McRAS-AC 
(Sud et al., 2013) used cloud water path scaling to account 
for the radiative effects of subgrid scale cloud water inho- 
mogeneity. More detailed discussions about RRTMG in the 
GEOS-5 AGCM can be found in Oreopoulos et al. (2012). 

2.3 Experimental design 

In order to investigate BBA effects on SEA climate, sev- 
eral observation-inspired experiments with and without BB 
emission over SEA are designed. Figure la shows the cli- 
matological amount of carbonaceous aerosol emission from 
BB during FMA averaged over 8 years from 2003 to 2010, 
and Fig. lb depicts the time series of the boxed area. The 
QFED data set described earlier was used for this figure. 
Massive BB emissions occur during FMA in the eastern re- 
gions of Cambodia, Myanmar, Laos, and northern Thailand 
with peaks in March. Figure lb shows large temporal vari- 
ations of BB emissions. It could be high in one month and 
also could become near zero in another month. 

To isolate the potential BBA effects, the AGCM exper- 
iments were performed with climatological SSTs to elimi- 
nate large-scale forcing (e.g., El Nino events) influences due 
to SST variability. Moreover, to separate the signal from 
the model’s own internal variability, multi-member ensemble 
simulations were performed. Each simulation-set consists of 
a 10-member ensemble covering the early January to late Au- 
gust period; each runs starts with different initial conditions 
taken from the model runs used in Sud et al. (2013). To es- 
timate the signal of aerosol effects on climate variability, we 
conducted experiments with “Zero” BB emission over the 
green dash box region (Fig. la) and compared them to exper- 
iments with “High” BB emission in 2007. BB emissions out- 
side of the boxed area and all other sources of aerosol were 
set to climatological means for both “High” and “Zero” emis- 


sion experiments. Such an “aerosol on” and “aerosol off’ de- 
sign is often used in aerosol-climate sensitivity studies (e.g., 
Lau et al., 2006; Wu et al., 2013) to better depict the aerosol 
signal, but this methodology has the drawback of making the 
comparison of simulations with observations difficult. Initial 
simulation sets using year-to-year emission data set did not 
yield statistically significant differences on circulation, while 
some sensitivity to enhanced emissions could be discerned 
in increased AODs, brightened liquid clouds and decreased 
rainfall. Clearly, by design, the differences between “High” 
and “Zero” emission experiments yield the effect of BBA, 
black carbon, organic carbon, and sulfate originating from 
the boxed area and in order to isolate statistically significant 
signals a Student’s t test is employed. We considered the dif- 
ferences exceeding the 95 % confidence level in a difference 
field as statistically significant. 

The differences between “HighBoth” and “ZeroBoth” 
simulations are a measure of the total BB effect. Here, 
“Both” means that the model’s experimental setup includes 
both aerosol direct and indirect effects. The indirect-only 
simulations are denoted by “Highlnd” and “Zerolnd” exper- 
iments. In these simulations, we neglect the aerosol direct 
effect by turning off aerosol radiative interactions globally, 
thereby allowing only the indirect effects of aerosols to op- 
erate on clouds and influence their radiative effects. Highlnd 
minus Zerolnd differences therefore measure the strength of 
the BBA indirect effect. Finally, while not additive due to 
nonlinearities, comparing “Both” and “Ind” runs gives in- 
sight into the relative contributions of direct and indirect ef- 
fects to the total aerosol effect. 
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3 Results 

3.1 Comparison between model simulations and 
satellite observations 

Some insight on the effects of BB in SEA can be obtained 
by comparing high and climatological BB conditions in ob- 
servations. We have constructed composites of MOD1S Aqua 
Level 3 daily products (MYD03D8) for 36 days of the high- 
est BB emission index within the 2003 to 2010 period. The 
BB emission index is defined as the area average for the 
dashed box area in Fig la. Smoothed time series of the index 
is shown in Fig. lb. Figure 2 shows comparisons between the 
high-emission days and the 8-year climatology of MODIS- 
Aqua AOD (Fig. 2a) and liquid cloud effective radius ( /Air, 
Fig. 2b). When BB is high in FMA over inland areas of SEA 
compared to normal days, anomalously high AOD appears 
over the northern part of SEA up to the coast of southern 
China. According to Lau and Kim (2014), low-level wind 
in the area carries BBA from the source region to southern 
China, resulting in the high BB AOD anomaly of Fig. 2a. Ac- 
cording to the rightmost panel of Fig. 2b, the corresponding 
negative anomaly of R e ff coincides with the region where a 
positive AOD anomaly exists. This is a classic manifestation 
of aerosol indirect effect whereby increased BBA reduces 
the size of cloud droplets by increasing CCN number con- 
centration. Indeed, if the negative anomaly of i? e ff is related 
to aerosol, then the imprints of other aerosol indirect effects 
may also exist in other meteorological fields, such as the pre- 
cipitation. Figure 2c compares composited daily GPCP pre- 
cipitation for the enhanced BB days and the climatology. The 
difference plot reveals a negative precipitation anomaly over 
the aerosol source and its adjacent areas while in the vicinity 
of the east coast of China increased rainfall is observed. This 
can be interpreted in a Lagrangian framework, by the cloud 
holding more cloud water due to reduced autoconversion ef- 
ficiency, but as the cloud advects downwind (i.e., towards the 
northeast direction), it eventually releases cloud water as pre- 
cipitation far away from the source region. Fligher BB emis- 
sion days are selected from every year and compared with 
climatology to remove interannual variability and SST forc- 
ing and isolate the BBA effects. Flowever, since large-scale 
SST forcings, such as El Nino, can simultaneously trigger a 
reduction in precipitation and an increase in aerosols (Tosca 
et al., 2010), it is necessary to study shifts in the precipitation 
pattern from AGCM simulations. 

In order to evaluate the sensitivity of the AGCM to BBA 
variations, model output is compared to the satellite data 
analysis. Figure 3 shows the overall BBA-simulated sensi- 
tivity of the model as the difference between HighBoth and 
ZeroBoth experiments. Anomalies in /Air are obtained using 
COSP’s MODIS simulator in the GEOS-5 AGCM for fair 
comparison with observations. For the 10-member ensem- 
ble mean, simulated AOD increases downwind (i.e., towards 
the northeast direction) of the BB source and correspond- 



Figure 2. Composite analysis and differences (Diff) in (a) aerosol 
optical depth and (b) liquid cloud effective radius (pm) from 
MODIS-Aqua retrievals, and (c) GPCP precipitation (mm day -1 ). 
High emission (left), climatology (middle) and differences: high mi- 
nus climatology (right) panels respectively. Here 36 high-emission 
days and 8-year climatology data are used. 


ingly R e ff decreases. Not only are McRAS-AC AGCM sim- 
ulations capable of simulating the response of cloud droplets 
to AOD, but the model’s overall response exhibits reason- 
able spatial coherence with the composite maps of observa- 
tion data (Fig. 2), particularly on the downwind side of the 
BB source. Wind vectors at 800 hPa are plotted on Fig. 3a 
to explain an advection of BBA. For precipitation, the obser- 
vations show a dipole-like anomaly pattern, namely overall 
decrease near source area and at eastern locations, and in- 
crease further east near the coast (Fig. 2c, Diff). In the model, 
on the other hand, the average BB signal on FMA precipita- 
tion materializes as a reduction in the precipitation of large 
areas south and east of the source region. Since SSTs were 
prescribed climatologically, less meaningful responses over 
the ocean is expected. Potential BBA effects on SEA pre- 
monsoon are investigated further in the following sections by 
taking appropriate differences of GEOS-5 experimental sets 
described earlier. 

3.2 BB effects on cloud microphysics and precipitation 
simulation 

One of the mechanisms that changes the nature and amount 
of precipitation is cloud microphysical processes as influ- 
enced by aerosols, widely known as the aerosol second in- 
direct effect (Albrecht, 1989). This mechanism acts on the 
autoconversion rate that modulates the intensity of liquid pre- 
cipitation. Our area of focus where precipitation decreases 
in MODIS analysis and model simulations is likely affected 
by BBA that are transported to southern China, where a per- 
sistent cloud band exists. In the model simulation the hori- 
zontal and vertical location of aerosol and the cloud band(s) 
are in close proximity, as seen in Fig. 4. The figure shows 
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Figure 3. Same variables as in Fig. 2 but for HighBoth minus ZeroBoth differences of AGCM simulations. Panels show (a) aerosol optical 
depth, (b) liquid cloud effective radius (|im) and precipitation rate (mm day - 1 ) during the FMA time period. Vectors on (a) are wind vectors 
at 800 hPa from HighBoth, only plotted where AOD increases. 
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Figure 4. Vertical cross section of simulated mixing ratio of BBA 
(shading, pgkg -1 ) and cloud liquid water (contour, mg kg -1 ) in 
March in HighBoth simulations; zonal averaging is performed for 
105 to 120° E. 


the vertical cross section of the BBA mixing ratio (shading) 
and cloud liquid water content (contour) for March, obtained 
from the HighBoth experiment in the vicinity of decreased 
precipitation as shown in Fig. 3c (105 to 120° E). BBA are 
lifted aloft by topography (oriented in a north-south direc- 
tion), and act as an additional source of CCN in pre-existing 
clouds which are mostly low-level (warm) and therefore of 
liquid phase at this particular location and time of the year. 

In the model, BB produces sulfate and carbonaceous 
aerosols that can be activated as cloud droplets. Sulfate 
aerosols are highly soluble, while the carbonaceous aerosols 
have both hydrophobic and hydrophilic modes. Hydrophilic 
organic and black carbon aerosols have prescribed fractions 
of soluble mass (0.25 and 0.1) so they also act as CCNs. 
Thus BBAs present at the level of developing liquid clouds 
also activate along with the background aerosols. Under con- 
ditions of massive BB, the CCN number concentration in- 
creases greatly, yielding increased cloud drop number con- 
centration and reduced cloud droplet sizes for constant cloud 
water. The underlying physics leading to the reduction in R e ff 
is well captured by the model (Fig. 3b). Smaller droplets re- 
duce the efficiency of autoconversion from cloud liquid wa- 
ter to rain, resulting in less precipitation. The double mo- 
ment microphysics in McRAS-AC reduces autoconversion 


rate process via the parameterization. 



auto 


-kl 4 c n~ 2 , 


( 1 ) 


where L c is the cloud liquid water content (kgm -3 ), N c is 
the cloud drop number concentration (m -3 ), and K is an ac- 
cumulated constant for autoconversion (see Eq. A. 2 in Sud 
and Lee, 2007, with units of kg -3 m 3 s -1 ). From Eq. (1), as 
N c increases under the assumption of constant L c , the au- 
toconversion rate decreases. The occurrence of this second 
aerosol indirect effect is shown by the model experiments. 
Monthly mean difference fields between HighBoth and Zer- 
oBoth experiments from March to May are shown in Fig. 5 
for AOD, N c , cloud liquid water path (LWP), precipitation 
and total cloud fraction. N c has been vertically averaged from 
900 to 750 hPa. Red (blue) color indicates positive (negative) 
anomaly with increasing BBA. Green contours delineate the 
areas where the change by BBA is significant at the 95 % sig- 
nificance level, based on Student’s t test. Aerosols clearly in- 
crease due to BB emission with an annual peak in March, and 
so does the AOD anomaly. Since February is a dry season for 
the area, the analysis focuses on March, April and May, with 
the latter month delineating the onset of the East Asian mon- 
soon. With BB occurring mainly in early spring, aerosol con- 
centrations in May should not be affected much, so any sig- 
nal in the meteorological fields for that month is potentially 
due to circulation and land surface changes induced by BB 
emissions in the preceding months. In the model we can de- 
compose AOD anomaly by species. The difference of AOD 
between HighBoth and ZeroBoth experiments over the boxed 
area in Fig. 1 for the month of March is 0.6, 0.487 com- 
ing from organic carbon, 0.063 from sulfate, and 0.05 from 
black carbon. There is some amount of “background” sul- 
fate in ZeroBoth (AOD = 0.129), but organic and black car- 
bon aerosols are mostly from BB emissions, so their AODs 
in ZeroBoth are quite small, 0.031 and 0.013, respectively. 
Dust and sea salt aerosol presence are very small over the 
region with AODs less than 0.01. 

As stated earlier, as the BBA loadings increase, grid mean 
N c , the product of in-cloud N c and cloud fraction also in- 
creases. Overall, both N c and LWP increase in the high 
BB experiments due to delayed precipitation, particularly in 
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Figure 5. HighBoth minus ZeroBoth representing BBA effects on 
(a) aerosol optical depth, (b) grid mean cloud drop number con- 
centration, (c) liquid water path (LWP), (d) precipitation and (e) 
total cloud fraction (%) from COSP. Green contour mark regions of 
>95 % significant in a Student’s t test. 


March and April. Still, some regions exhibit negative grid 
mean N c and LWP anomalies, possibly because of reduced 
cloud fraction (Fig. 5e) due to reduced grid-scale relative 
humidity (RH) that detennines cloud amount for stratiform 
clouds. The reduced RH is the outcome of larger stability of 
the lower atmosphere, which suppress rising motion. 

While enhanced BB emission increases AOD, CCN, N c , 
and even L c , the relationship is not linear. Increased L c 
eventually creates a tendency for higher autoconversion and 
precipitation rates, which opposes the tendency of the in- 
creased N c (Eq. 1). If we do not account for complex feed- 
backs, precipitation near the BB emission source can be 
expected to decrease if the increased N c effect is stronger 
than the enhanced L c effeet, as shown in Fig. 5d for March 
and April. While the satellite data analysis suggests alternat- 
ing negative-positive precipitation anomalies along the wind 
flow, a weak positive anomaly surrounds the simulated strong 
negative anomaly over the South China Sea and northern 
China in April (Fig. 5d). This may be explained by two 
possible mechanisms. Liquid cloud water gets transported 
downstream instead of precipitating out locally because of 
suppressed autoconversion, while reduced local precipita- 
tion creates favorable circulation conditions for precipitation 
downwind. Meanwhile, a statistically significant anomaly of 


precipitation is found in May, suggesting that large BB emis- 
sion in March and April can have a delayed effect even in 
regions far away from the source. Microphysical processes 
may therefore not be the only mechanism that reduces pre- 
cipitation. The impact of BBA on May precipitation could be 
a combination of direct, indirect effects, and feedback pro- 
cesses initiated by aerosols in March and April. Further anal- 
ysis of circulation changes is needed to distinguish whether 
this anomaly can be attributed to cloud microphysics or some 
other mechanism, a topic that we will address in Sect. 3.4. 

3.3 BB effects on the radiation budget 

BBA can change the radiation balance by both their direct 
and indirect effects. The direct effect of BBA consists of 
scattering (sulfate aerosols) and absorption (black carbon 
aerosols) of incoming solar radiation which cause surface 
cooling and atmospheric heating. As discussed in Sect. 3.2, 
the indirect effect of BBA comes from altering cloud optical 
properties like A* c ff and cloud amount which modify the net 
(= shortwave + longwave) radiation budget at the top of the 
atmosphere (TOA), atmosphere (ATM), and surface (SFC). 
Figure 6 illustrates the magnitude of the net radiation change 
at TOA, ATM, and SFC due to both direct and indirect 
aerosol effects, primarily due to changes in shortwave (SW) 
radiation. Each map shows the monthly mean difference be- 
tween HighBoth and ZeroBoth experiments from March to 
May with red (blue) indicating heating (cooling) anomalies 
by aerosol, and green contours delineating the areas of statis- 
tically significant change. The overall net radiative effect of 
BBA is TOA/SFC cooling, and ATM heating near the source 
region, but its interpretation requires further scrutiny because 
contributions to net radiation change also come from circula- 
tion changes and associated feedbacks. The ATM heating in 
March and April provides a clearer signal of direct effects 
since the radiative heating comes almost exclusively from 
aerosol absorption, while TOA and SFC cooling comes from 
both direct and indirect aerosol effects. In May there is little 
aerosol direct effect, but a significant anomaly signal exists at 
the TOA and SFC due to feedbacks from circulation changes. 
When examining TOA and SFC radiation fields in May, the 
dipole pattern seen near the east coast of China and Korea 
is due to cloud fraction change (Fig. 5e) consistent with the 
precipitation change shown in Fig. 5d, and discussed further 
in the following section. 

Table 2 shows the differences between HighBoth and Ze- 
roBoth experiments of the net downward (down minus up) 
fluxes in March when BBA peaks, regionally averaged across 
the emission control region. Aerosol radiative effects are 
much larger in the SW than the longwave (LW), so most of 
the net radiation change comes from SW effects. The cor- 
responding clear sky fluxes demonstrate that BBA increases 
SW reflectance, but also absorbance, because large fractions 
of BBA are composed by carbonaceous aerosols, which are 
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Figure 6. Same layout as in Fig. 5, but for net radiative fluxes at (a) 
top of the atmosphere, (b) column atmosphere, and (c) surface. 

Table 2. Radiative flux change (W m -2 ) by BBA: numbers indi- 
cate the difference between HighBoth and ZeroBoth experiments in 
March and April, regionally averaged from 90 to 1 10° E and from 
12 to 30° N only including land grid (emission control region). All 
fluxes are net downward, which means upward fluxes are subtracted 
from downward fluxes. 



TOA 

ATM 

SFC 

SW, all sky 

-9.5 

15.1 

-24.6 

SW, clear sky 

-9.0 

17.1 

-26.1 

LW, all sky 

0.3 

-2.5 

2.8 

LW, clear sky 

1.2 

-1.0 

2.2 


efficient absorbers of SW radiation. CRE is defined as 

CRE = /'all— sky E c i ear _ s ky, (2) 

where F is the net downward flux at the TOA or surface. 
The SW CRE change by BBA (the difference between all 
sky and clear sky in Table 2 - not shown) is, somewhat sur- 
prisingly, positive at the surface (weaker SW CRE for the 
HighBoth experiment). Despite LWP increases (Fig. 5c) and 
R e ff decreases (Fig. 3b) in conditions of high BB, yielding 
an average optical thickness increase of 46 % for the cloudy 
part of the BBA source region, decreased total cloud frac- 
tion (Fig. 5e) due to circulation changes overcomes increased 
cloud brightness. So, the total indirect effect, namely R e ff de- 
creases (classic “Twomey effect”), accompanying LWP in- 
creases and any resulting cloud feedbacks, seem to counter- 
act the direct aerosol effect in these GEOS-5 experiments. 

3.4 Temperature, moisture, and circulation changes 

In this section we discuss meteorological the consequences 
of the BBA radiative effects, which we have shown to lead to 
surface cooling and atmospheric heating. The surface tem- 


perature anomaly due to BBA is plotted in Fig. 7 and can be 
seen to be very significant near the source region in March 
and April. Weak negative anomalies also appear in May, 
but are mostly insignificant statistically. The vertical profile 
of temperature change (black line) by BB, regionally aver- 
aged from 100 to 120° E and from 18 to 30° N (cf. red box 
in Fig. 7), is plotted in Fig. 8. This profile is obtained as 
the difference between the HighBoth and ZeroBoth exper- 
iments in March and April when the decrease of precipi- 
tation is significant, and reveals the presence of a cooling 
signal from the surface all the way up to the 250 hPa level. 
In order to better understand what causes the temperature 
change, the model’s major heating/cooling rate contributions 
are shown in Fig. 8. The orange line shows the SW heat- 
ing rate (Kday -1 ) anomaly, the red line the LW heating rate 
anomaly, and the blue line the anomaly of the heating rate 
due to the model’s moist physics, namely large-scale con- 
densation and convective processes. As expected, SW radi- 
ation heats the atmosphere near the height of the aerosol 
layer (Fig. 4). The reason the temperature profile does not 
cross over to the positive side is due to other contributors to 
temperature change, namely LW and moist physics, both of 
which cool the low and middle troposphere. The increase in 
LW cooling is a consequence of increased cloud liquid water 
between 800 to 600 hPa due to aerosol-induced changes in 
cloud microphysical processes. Although the magnitude of 
LW cooling is only about a quarter of the SW heating, it has 
an impact since LW cooling occurs at the time and location 
of SW heating, because of liquid water and BBA collocation. 

The major factor contributing to the negative temperature 
anomaly is the reduced moist heating, the most significant 
change of all the heating rate components of the model’s 
physics. A negative moist physics heating rate anomaly 
translates to subdued cloud formation by large-scale con- 
densation and even moist convection. In the area of interest, 
March precipitation mostly comes from large-scale conden- 
sation, while in April there is some contribution from moist 
convection. The reduced convective precipitation in April, 
accounting for about 40 % of the total precipitation reduc- 
tion, can be explained by changes in the vertical tempera- 
ture gradient. BBA direct radiative effects make the surface 
cooler and the 700 hPa level warmer, thus decreasing low- 
level atmospheric instability as seen in the vertical temper- 
ature profile anomaly; even though the overall temperature 
change is negative (cooling), a bump of temperature anomaly 
forms near 700 hPa that suppresses the onset of moist con- 
vection. 

Another important reason behind moist physics suppres- 
sion, which can explain both large-scale condensation and 
moist convection, is change in atmospheric moisture con- 
tent. Zonally averaged moisture and meridional circulation 
anomalies due to BBA within 100-120° E for March and 
April, and 1 10-140° E for May are plotted in Fig. 9. The 
blue shading, indicating dry anomaly, spreads over the 20- 
30°N latitude zone where the BBA sources are located. A 
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Figure 7. Same layout as in Fig. 5, but for surface temperature (K). Red dashed domain is for area-averaged fields in Fig. 8 and Table 4. 



Figure 8. Vertical profile of temperature (K) and heating rate 
(K day - 1 ) differences between HighBoth and ZeroBoth simulations 
for March and April when precipitation decrease is significant; the 
averaging region is 100 to 120° E and 18 to 30° N and is marked 
in Fig. 7. Black, orange, red and blue lines represent temperature, 
SW heating rate, LW heating rate, and the heating rate due to the 
model’s moist physics, respectively. 


few factors play a role in making this region drier. One is 
reduced surface evaporation (Fig. 10) in the region of neg- 
ative surface temperature anomaly. The other is circulation 
changes, specifically the substantial downward and south- 
ward flow anomalies induced by BBA. While the down- 
ward anomaly could be the result of reduced moist activ- 
ity, the accompanying southward anomaly may actually be 
the cause of reduced moisture transport from low latitudes. 
The column-integrated moisture convergence anomaly (not 
shown) in the region where precipitation decreases is nega- 
tive with some degree of statistical significance, albeit less 
than surface evaporation (Fig. 10). Another possible cause 
for overall drying is the decreased precipitation itself, imply- 
ing positive feedback. Because of the BBA indirect effect, 


reduced autoconversion leaves behind more in-cloud water, 
which advects downwind instead of being converted into lo- 
cal precipitation, resulting in a reduced supply of moisture 
to the levels underneath and creating a feedback loop where 
weaker latent heat flux at the surface causes further decreases 
in precipitation. 

Several mechanisms can potentially reduce precipitation 
in the downwind side of an active BB region. Cloud micro- 
physics can delay the precipitation process by slowing down 
autoconversion, and then radiation can help make the area 
stable and dry, all conditions unfavorable for vigorous moist 
processes. Moreover, dry anomalies can be the result of cloud 
microphysics as well as in lack of rain and its evaporation. 
Likewise, a number of other variables may be changing in 
the same direction due to direct BBA effects on radiation 
and the indirect effects on cloud microphysics. Both effects 
cause SW dimming at the surface, low-level drying, and de- 
creased precipitation. Separating microphysical from radia- 
tive effects is thus a worthwhile objective which we pursue 
in the following section. 

3.5 Quantitative breakdowns of direct and indirect 
effects 

In the previous sections, all the results explaining aerosol ef- 
fects were based on FlighBoth and ZeroBoth experiments, 
the first including BBA from a high-emission year and the 
latter entirely neglecting BBA emissions from the area of 
strongest fire activity. The GEOS-5 AGCM accounted for 
both direct elfects in its radiative transfer routines and in- 
direct effects in its cloud microphysics routines. Differences 
between the two experiments capture both aerosol direct and 
indirect effects (as well as feedbacks), in other words, com- 
bined effects (CE). In two other experimental sets, the ‘High- 
Ind’ and Zerolnd experiments, direct effects of aerosol on 
radiation are ignored (globally), leaving only the indirect ef- 
fect (IE) of BBA to be diagnosed as the difference between 
the two Ind experiments. Our diagnostic approach to sepa- 
rate the direct and indirect aerosol effects of a rather complex 
regional climatic response consists of comparing key vari- 
ables from “CE” and “IE” differences, both including feed- 
back from circulation changes. 
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Figure 9. Zonally averaged profiles of moisture (shading) and meridional circulation anomalies (vectors, horizontal component is for merid- 
ional wind anomaly, vertical component is for pressure velocity) from HighBoth minus ZeroBoth experiments over the longitude sector 
100-120° E for March and April, and sector 1 10-140° E for May. Units of pressure velocity, meridional wind, and water vapor mixing ratio 
are 10 -2 Pas - , ms -1 , and g kg - 1 , respectively. 


Surface Evaporation change [mm/day] 

High minus Zero 



Figure 10. Similar to Fig. 5, but for surface evaporation (mm day 1 ). 


Table 3 shows the radiative fluxes in the same way as Ta- 
ble 2, but for “IE” and “CE minus IE”. Evidently, CE of TOA 
and SFC SW fluxes and atmospheric column SW absorption 
are much larger than the corresponding IE of aerosols on SW. 
This implies a much stronger contribution of the direct effect 
(DE) of BBA in CE and makes sense because BBAs have 
large AOD over the high-emission regions. Moreover, BBA 
effects on radiation in the CE runs are quite similar for clear- 
sky and all-sky conditions, as pointed out earlier in Sect. 3.3. 
This provides further support for the notion that the IE due 
to BBA is an order of magnitude smaller on the SW and net 
radiation compared to the corresponding DE of aerosols as a 
component of CE. Net radiation change by IE turns out small 
because it depends on both cloud fraction (which depends 
on cloud production, cloud dissipation, and cloud advective 
tendencies) and cloud optical thickness (which depends on 
CCN and cloud water removal by precipitation). In Zerolnd, 
the (BBA-independent) cloud fraction increases while the 
cloud optical thickness decreases compared to Highlnd sim- 
ulations. 

One of the interesting features of the BBA signal is de- 
creased precipitation over the downwind side of the source. 
The separation of “CE” and “IE” impacts on precipitation 
would be interesting to study for this area. To minimize 
feedback contributions, we focus on variables that are pri- 
marily forced directly during the March and April time 
frame and near the source region, in particular 100 to 120° E 
and 18 to 30° N. Table 4 provides the spatiotemporal aver- 


Table 3. Radiative flux change (W m -2 ) by indirect effect of BBA, 
and the differences from Table 2. “IE” indicates the difference be- 
tween Highlnd and Zerolnd experiments, “CE” indicates the differ- 
ence between HighBoth and ZeroBoth experiments, while the oth- 
ers are the same as Table 2. 


TOA ATM SFC 



IE 

CE-IE 

IE 

CE-IE 

IE 

CE-IE 

SW, all sky 

0.5 

-10 

0.0 

15.1 

0.5 

-25.1 

SW, clear sky 

-0.3 

-8.7 

-0.3 

17.4 

0.0 

-26.1 

LW, all sky 

-2.1 

2.4 

0.8 

-3.3 

-1.3 

4.1 

LW, clear sky 

-0.4 

1.6 

0.7 

-1.7 

-1.1 

3.3 


ages of these CE and IE breakdowns. While the CE pre- 
cipitation reduction in HighBoth minus ZeroBoth BBA is 
1.08 mm day -1 , the corresponding IE precipitation reduction 
is 0.77 mm day -1 only. For a linear system, one would at- 
tribute the 0.31 mm day -1 reduction corresponding to the CE 
minus IE difference, to the direct aerosol effect, but we are 
well aware that linearity is not necessarily a good assump- 
tion, so we view the differences as representing add-on di- 
rect effects that also contain effects of interactive circulation 
changes. Even though the IE averages do not show much 
change in the simulated surface temperature and evapora- 
tion of the boxed region, IE does have a prominent role in 
decreasing surface precipitation, which is caused not only by 
autoconversion reduction, but also by low-level drying due to 
SW dimming in cloudy areas. In other words, the suppressed 
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Table 4. Analysis of combined effects (CE, Direct+Indirect) and 
only indirect effect (IE). 


Aerosol 

effect 

Combined 

effect 

Indirect 

effect 

Precipitation 3 (mm day -1 ) 

-1.08 

-0.77 

Surface temperature 3 (K) 

-0.59 

0.09 

Surface evaporation 3 (mm day - 1 ) 

-0.25 

0.07 

Moisture b (gkg -1 ) 

-0.38 

-0.23 

Moist heating rate b (K day - 1 ) 

-0.29 

-0.16 

Temperature 4 5 * * * * * 11 (K) 

-0.24 

-0.07 


a Regional average from 100 to 120° E and 18 to 30° N, on March and April, 
and b 975 to 500 hPa vertical. 


autoconversion that follows CCN and N c increases due to 
the presence of additional BBA in the IE simulations and 
decreases precipitation, which creates a dry anomaly in the 
atmosphere beneath the precipitating cloud due to reduced 
evaporation of rain. In comparison, the DE as a part of CE 
has a more straightforward effect on the moisture supply that 
can be traced to atmospheric stabilization and reduced sur- 
face temperature due to surface cooling. So while both CE 
and IE tend to reduce precipitation, the mechanisms can dif- 
fer overall, despite sharing the common processes of slower 
autoconversion and low-level drying. 


4 Summary and discussion 

An aerosol impact study including both the direct and in- 
direct effects focusing on the Southeast Asia pre-monsoon 
season is conducted based on simulations using the GEOS- 

5 AGCM with double moment cloud microphysics called 
McRAS-AC, interactive GOCART aerosol model, advanced 
radiative transfer package RRTMG applying the Monte Carlo 
Independent Column Approximation, and the CFMIP Obser- 
vation Simulator Package (COSP). Analysis of GEOS-5 in- 
tegrations with and without BB emission allows us to sepa- 
rate the responses of clouds and precipitation to aerosol from 

those due to changes in meteorological fields. Our analysis 

indicates that plausible reasons for the reduced precipitation 

are (a) vertical stabilization by atmospheric heating aloft ac- 
companied by surface cooling due to the shortwave scattering 

and absorption by the BBA; (b) less efficient autoconversion 
despite liquid water increases due to increased cloud droplet 
number concentration; and (c) suppressed moist processes 

due to atmospheric drying. With properly designed experi- 
ments we managed to separate the impacts of direct and in- 

direct effects. While vertical stabilization is traced to direct 
aerosol-radiation interaction, which causes rapid cloud ad- 
justments (commonly referred to as the “semi-direct effect”) 
because of depressed dynamical forcing, and the reduced 
autoconversion rate is primarily a consequence of aerosol- 


a) Temperature b) Zonal wind 



Figure 11. Zonal mean temperature (K) and wind (ms -1 ) differ- 
ences between HighBoth and Zero Both BBA for 1 10-140° E in 
May. Green contour identifies regions with > 95 % significant dif- 
ferences according to the Student’s t test. Black contour on (b) is 
zonal mean wind in HighBoth run. 


cloud interaction (the indirect effect), the drying of the lower 
and middle troposphere is caused by both. 

An interesting, and somewhat unexpected, consequence of 
enhanced BBAs is the May precipitation anomaly near the 
Korean peninsula shown in Fig. 5d. Since BB is not a ma- 
jor factor in May aerosol loadings, the precipitation anomaly 
could be due to circulation and land surface changes trig- 
gered by BB in the preceding month. In March and April, the 
surface temperature over Southeast Asia drops significantly 
due to the combined direct and indirect solar dimming effect 
of BBA and this reduces the meridional temperature gradient. 
Figure 9 shows that the overall circulation anomaly heads 
south in March and April. The May circulation anomaly ex- 
hibits downward motion at 30° N and a little upward motion 
south of 30° N. According to Kim et al. (2007) an upper level 
jet stream change can induce secondary circulation changes 
near the entrance of the jet core in East Asia. In their analysis, 
an initial surface cooling by the direct effect of sulfate aerosol 
results in a reduced north-south thermal gradient. Figure 11a 
shows a similar weakened meridional temperature gradient 
change due to surface cooling found in March and April. 
This reduced gradient weakens the zonal wind shear through 
the thermal wind relationship, and slows down the westerly 
jet stream (Fig. lib). The deceleration causes ageostrophic 
meridional winds and, in this case, anomalous sinking mo- 
tion at 30° N (see Fig. 9, May), conditions that are less favor- 
able for precipitation. Although Kim et al. (2007) account 
only for direct forcing of aerosol, the circulation anomalies 
induced by the BBA emissions of this study are similar, be- 
cause indirect effects did not affect much the surface forcing 
in CE simulations. 

While this study provided some confirmation that our BB 
sensitivity in the model looks similar to that from the MODIS 
analysis, the “Zero” BB assumption is admittedly extreme. 
So the year-to-year change of meteorological fields by BBA 
could be weaker than suggested by the results shown here. 
But given the plausibility of how the model’s mechanisms 
operate, there is good possibility that real conditions would 
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be consistent with the overall tendencies of the model. Still, 
there is much room for further development of the GEOS-5 
model towards more realism. For example, phenomena such 
as aerosol-induced convective invigoration (Rosenfeld et al., 
2008) cannot be properly reproduced in our model because 
heat release due to freezing does not affect the convective 
mass flux. This process could be better represented in a bulk 
mass flux convection scheme (e.g., Kim and Kang 2012), but 
it remains to be seen whether its inclusion in such a scheme 
would ultimately affect overall convective activity in a sub- 
stantial way. Evidently, further refinements and cloud model 
validations are needed for a better understanding of the role 
of aerosol-convection interactions in the seasonal develop- 
ment of the summer monsoon. Our method of separating di- 
rect and indirect aerosol effects may be imperfect, but no bet- 
ter alternative currently exists given present modeling limita- 
tions. Regardless, we believe that this study provides a foun- 
dation on which to develop better methodologies to properly 
distinguish direct and indirect effect sensitivity to aerosols in 
large-scale models. 
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